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Abstract. Detailed models are presented for the late epoch mid infrared (MIR) emission from collisionally heated 
grains in the shocked circumstellar gas around SN 1987A. Thermal dust emission from a region of moderate density 
interior to the thick inner ring seen with the Hubble Space Telescope (HST) is found to be a natural explanation 
for the MIR spectral energy distribution measured by ISOCAM. The MIR-spectrum can be reproduced by a 
mixture of silicate-iron or silicate-graphite grains or by a composition of pure graphite grains. A composition of 
pure iron grains on the other hand can be excluded and a pure silicate composition does not seem to be very likely. 
The dust-to-gas ratio in the interaction zone is ~ 0.01%, an order of magnitude lower than estimates for dust 
abundances in the winds of red supergiant (RSG) stars in the LMC. This low dust abundance can be accounted 
for by a combination of evaporation through the UV-flash from the supernova outburst and subsequent sputtering 
in the shocked gas. For this explanation to hold, dust in the pre-supernova circumstellar medium (CSM) would 
have to have been predominantly composed of grains other than graphite, with a maximum size smaller than 
~ 0.1 /im. 
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1. Introduction 

SN 1987A, a supernova of Type II, made it possible for 
the first time to make detailed observations of the inter- 
action of a supernova ejecta with the very innermost re- 
gion of its CSM. IR measurements enable us to analyse 
the abundance, composition, and size distribution of the 
circumstellar dust grains and to study grain destruction 
processes such as evaporation by the UV-flash of the su- 
pernova outburst and sputtering in the shocked gas behind 
a very strong shock. From this, new insights into the his- 
tory of the progenitor star can be obtained, providing a 
more complete picture of the SN 1987A. 

SN 1987A is the first supernova for which the progen- 
itor was observed prior to outburst. It has been identi- 
fied as the most luminous star of the system Sanduleak 
-69°202 in the Large Magellanic Cloud (LMC) (West, 
I1987[) which had been classified as a blue supergiant (BSG) 
of spectral Type B3 I (Rousseau. ll978fl . Observations with 
the HST showed a complex axis symmetrical structure of 
its CSM with a thick inner ring with a diameter of ~ 1'.'6 
and two outer rings with larger diameter at each side seen 
with a viewing angle of ~ 43° (Burrows et al. . 11995) ). The 
analysis of dust scattered light of the supernova outburst 
suggest that the rings are connected with gas and dust 



Send offprint requests to: Jorg Fischera 
Correspondence to: Joerg.Fischera@mpi-hd.mpg.de 



distributed in an hour glass - shaped shell (Crotts et al., 

mag. 

Initially, the supernova ejecta freely expanded in the 
thin wind of the BSG, driving a blast wave into the 
CSM with a velocity of ~ 30 000 km/s. The reappear- 
ance of the radio emission (Staveley-Smith et al.,[1992) 
almost coincidentally with the appearance of soft X-ray 
emission detected by ROSAT (Beuermann et al., 
Gorenstein et al.. I1994|l indicated that after ~ 1200 days 
the blast wave had reached denser material interior to the 
thick inner ring, which slowed down the shock velocity 
to 2900 ± 480 km/s (Gaensler et al.,E&3- Chevalier & 
Dwarkadas l)1995f) supposed that this denser gas is com- 
prised of material from a RSG phase of the progenitor 
of SN 1987A before it evolved into the BSG that finally 
exploded. They refered to this region, in which the gas 
should be ionised through the photon flux of the progeni- 
tor star of SN 1987A, as the "HH-region" . The interaction 
between the fast thin wind of the BSG and the slow mov- 
ing thick wind of a RSG is also thought to be responsible 
for the larger structure of the CSM (Blondin & Lundqvist, 
IT9*9*51 Martin & Arnett. HMSjl . 

In a previous paper (Fischera et al., [2002, paper I) 
we presented MIR measurements, made with ISOCAM 
(Cesarsky et al. 119960 towards SN 1987A 11 years after 
the outburst. These reveal the central region around the 
supernova position as a resolved MIR source with an ex- 
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tension and orientation consistent with the elliptical pro- 
jection of the thick inner ring, suggesting that the MIR 
emission is mainly circumstellar in origin. We found this 
emission is most probably from dust, collisionally heated 
in the shocked gas downstream of the blast wave as it 
expanded into the material of the Hll-region interior to 
the thick inner ring. On a purely energetical basis, all the 
emission could be from condensates in the metal rich core 
region of the expanding ejecta fFischera. 12000(1 . However, 
as argued in paper I, an emission which arises mainly from 
condensates does not explain the measured extension or 
orientation of the MIR source. 

In this paper we present detailed calculations of emis- 
sion from grains in the shocked CSM to analyse the im- 
plications of the ISOCAM measurements for the abun- 
dance, composition and size distribution of the circum- 
stellar grains. In Sect. we describe the dust model 
used to analyse the dust properties in the shocked CSM. 
Quantitative results are presented in Sect. 3. In Sect. 4 
we discuss probably the most important destruction pro- 
cesses, which are evaporation of grains through the UV- 
flash and sputtering in the shocked gas downstream of the 
blast wave. In Sect. 5 this is used to derive information 
on grain abundance, and composition in the CSM prior 
to the supernova outburst. A summary is given in Sect. 6. 
As in paper I we will assume a distance of 51 kpc to the 
supernova. 

2. Modelling the MIR emission from the shocked 
CSM 

We consider grains that are collisionally heated in the 
shocked circumstellar plasma downstream of the blast 
wave. 




Fig. 1. Cartoon of the basic structure caused by the in- 
teraction of the expanding ejecta with the HII-rcgion. 

The structure caused by the interaction of the expand- 
ing ejecta with its CSM is visualised in figure ^ The 
shocked circumstellar gas is compressed to a thin layer 
that grows approximately as O.lugi (Chevalier, 1982), 
where vs is the assumed constant speed of the blast wave 
and t the time since the shock reached the Hll-region. 



Downstream it is bounded by the contact discontinuity 
to the shocked outer parts of the expanding ejecta. This 
ejected material is heated by both the reverse shock and 
by a shock reflected from the inner boundary of the HII- 
region (Borkowski et al., 1997). It is thought to be the 
origin of the strong emission in Lya and Ha at a distance 
of ~ 0.6 of the radius of the thick inner ring (Michael et al., 
1998J, that was measured by the HST close to the time 
of the ISOCAM observations. We do not expect strong 
MIR emission from the shocked ejected material for the 
following reasons: 

1. The material at the reverse shock is mainly hydrogen 
and has such a low abundance in metals (e.g. Woosley, 
1988), that, if any, only very few grains may have been 
formed in this region of the ejecta. 

2. Further, as we show later (see Sect.EJ), the grain charge 
and the strength of the magnetic field in the shocked 
CSM in the interaction zone makes it unlikely, that an 
observable amount of circumstellar grains have been 
overtaken by the outer parts of the ejecta. 

3. The gas density of the ejecta at the reverse shock is, 
following the numerical results from Borkowski et al. 
I|1997[) . much less than the density downstream of the 
blast wave. In consequence, the heating of any ejected 
grains at the reverse shock should be much weaker than 
for circumstellar grains downstram of the blast wave. 

2.1. Gas parameters 

The collisional heating of a grain in a hot plasma depends 
mainly on the temperatures of electrons and ions and the 
number density of the different species of the gas. 1 We 
fix these properties to be broadly consistent with the ob- 
served radio and X-ray emission from the CSM. An analy- 
sis of the X-ray observations made with the ROSAT satel- 
lite towards SN 1987A up to and including the epoch of the 
ISOCAM observations is given in appendix [S] For tem- 
peratures above ~ 10 7 K, as encountered in the shocked 
gas, the grain heating of the smaller grains principally con- 
tributing to the emission observed by ISOCAM is mainly 
determined by the plasma density (see e. g. Dwek & 
Arendt 1992). As the physical conditions are somewhat 
uncertain we will consider two cases (model I and II) with 
different temperatures and densities. The parameters of 
model I and II are summarised in Table ^ and will be 
justified in the following. 

Borkowski et al. (1997) have shown that an assump- 
tion of constant density, as suggested by Chevalier & 
Dwarkadas, can explain the soft X-ray emission (Hasinger 
et al., 1996 ) until at least 3000 days after outburst. 
They derived a pre-shock hydrogen density of 75 cm~ 3 
with a corresponding shock speed «s of the blast wave of 
4100 km/s. In their calculations the Hll-region was mod- 

1 In general the heating, in particular the heating by the ions, 
has also a dependency on the relative motion of the grains in 
comparison to the motion of the gas, which is neglected for 
simplicity. 
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Fig. 2. Temperature distributions of spherical silicate and graphite grains for various grain sizes in the shocked 
CSM, assuming a density of the shocked gas of nn — 300 cm' 3 (model I). The variation in size between two lines is 
dlog(a[/xm]) = 0.1. The temperature distributions, shown with thicker lines, belong to grains with a radius of 0.001 /im 
and 0.01 /im. The steps in the p(T)-curves are due to the approximation used for ion heating and appear roughly at 
temperatures where the thermal energy of the grains is equal to the threshold energy of the ions (see also Popescu et 

ai- onnnji . 



Table 1. Parameters of the shocked gas 





model I 


model II 


shock velocity vs 


4100 km/s 


2900 km/s 


electron temperature T c 


2 x 10 7 K 


2 x 10 7 K 


ion temperature 71 


6 x 10 s K 


3 x 10 8 K 


hydrogen density uh 


300 cm -3 


600 cm -3 


helium density nn c 


2.5 (n Hc )© 


2.5 (n Hc )© 


metallicity Z 


0.3 Z Q 


0.3 Z Q 



elled as a thick torus in which the inner ring seen by the 
HST is embedded. They also mentioned that for a dif- 
ferent structure of the Hll-region higher densities would 
be possible. Indeed, a shock speed of 2900 ± 480 km/s 
as derived from radio observations (Gaensler et al.. H997[) 
would require a hydrogen density of the Hll-region of 
about 150 cm -3 . It might also be possible, that the den- 
sity increased with time (see appendix EJ although the 
X-ray measurements are also consistent with a constant 
external density. 

For simplicity we assume that at the time of the 
ISOCAM observations the blast wave was still expanding 
into a homogeneous Hll-region. We will consider hydrogen 
densities of 75 cm -3 (model I) and 150 cm -3 (model II) 
for the pre-shock region. The density of the gas is assumed 
to be compressed by a factor of four downstream of the 
blast wave and to stay constant in the whole interaction 
zone. 

In converting the densities given by Borkowski et 
al. into hydrogen number density and ionic abundances 
needed for the calculation of dust heating (Sect. 12. 2f) we 
adopted the abundances of the thick inner ring given 
by (Lundqvist & Fransson, I199fi|) . On this basis we took 
Helium to be a factor of 2.5 more abundant than in the 



sun, and approximated the metallicity of all heavier ele- 
ments to be 0.3 solar. 

Because of the long equipartition time of the shocked 
circumstellar plasma the electron temperature should be 
much lower than the temperature for the ions and should 
increase with the distance to the outer shock (see e.g. 
Burrows et al. l2000f) . Thus, in calculating the grain heat- 
ing we use two different temperatures for electrons and 
ions. As for the density, the temperature of the shocked 
CSM is taken to be independent of position. 

The ion temperature Tj is taken to be (see e.g. Longair, 
1997) 



7-1 n ra H 



(7 + I) 2 k B 



(1) 



with 7 = 5/3 and an atomic weight of fi = 1.6 appropriate 
to the abundances in the CSM. 

For the electrons we choose a temperature close to 
those derived from X-ray observations and predicted by 
numerical calculations appropriate to SN 1987A. From X- 
ray observations made with the ROSAT-satellite Hasinger 
et al. I|199fi|l derived a temperature of approximately T e « 
1.2 x 10 7 K. Analysing the X-ray spectrum of SN 1987A, 
taken later with CHANDRA, using a shock model with a 
constant T c , a higher electron temperature of the order of 
~ 3.5 x 10 7 K was found (Burrows et al. 120001 Park et al., 
I2()()2jl . In our calculation we will adopt T c = 2x 10 7 K. 

2.2. Calculation of the dust emission 

The dust in the CSM is taken to be spherical for simplic- 
ity as there is no observational evidence for other grain 
shapes in the CSM of SN 1987A. Spherical grains are also 
generally assumed in the literature for grains produced 
in stellar winds (e.g. Gail & Sedlmayr, 1999). However, if 
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the grains are not too different from being spherical this 
assumption should not influence significantly our results 
about the properties of the grains in the CSM. For exam- 
ple, it has been found, that the temperature of spheroidal 
grains heated by the interstellar radiation field (1SRF) can 
only vary by more than 10% if the axial ratio exceeds a 
value of 2 (Voshchinnikov et al. . 11999^1 . 

The grains are heated by collisions with electrons and 
ions of the hot ionized shocked plasma. The gas species 
is assumed to transfer all (if it sticks) or only a part (if 
it is not stopped) of the kinetic energy into thermal en- 
ergy of the grain. The calculation of the energy deposi- 
tion of non-stopping particles is based on their stopping- 
distances (ranges) in solids 2 . For the heating by ions we 
considered in addition to hydrogen and helium also the 
next most abundant elements oxygen, nitrogen and car- 
bon. The emission from larger grains we derived from their 
equilibrium temperatures. For smaller grains, where the 
deposited energy is typically larger than the thermal en- 
ergy, we took their temperature fluctuations into account. 
The model for stochastic dust emission from a hot plasma 
is described in more detail in Popescu et al. ( 2000 ) . As typ- 
ical dust species we consider graphite and silicate grains. 
The temperature distributions of very small silicate and 
graphite grains in model I are shown in Fig. El where it 
can be seen, that small grains can be heated to very high 
temperatures. The smallest grains with 10 A radius will 
reach temperatures well above their evaporation tempera- 
ture. The grains contributing to the measured IR emission 
should therefore not be smaller than this size. 

Because iron is expected to form in cool stellar out- 
flows independent of the C/O ratio and is potentially one 
of the main condensates in the circumstellar environment 
of oxygen rich stars (see e.g. Gail & Sedlmayr, 119991 and 
references there), we also carried out calculations for pure 
iron grains. For electrons in iron grains we used the analyt- 
ical expression for the electron range in graphite derived 
by Dwek & Smith (1996) on basis of observational data, 
correcting for the different density. The optical properties 
of iron spheres we derived using MIE-theory (Bohren & 
Huffman. I1983[) . whereby we included the dependence of 
the dielectric function in the IR on size and temperature 
of the iron grains (Fischera, 2000). For the heat capacity 
of iron grains cooler than 298 K we used values tabulated 
in the American Institute of Physics Handbook (1972). 
For iron grains warmer than this we used an analytical 
expression for the heat capacity given by Chase ([1998 ) , 
assuming that iron is in the a-(5-phase. 



2 The stopping power is due to ionisation losses in case of 
electrons. In case of the ions at energies considered here the 
stopping power is partly due to interaction with the electrons 
and partly due to inelastic scattering with the nucleons of the 
target material. 



2.2.1. Grain temperatures in the shocked CSM 

The variation of the equilibrium temperature of spheri- 
cal silicate, graphite, and iron grains in the shocked CSM 
is shown in Fig. [21 The different behaviour of the three 
species allows some information about dust composition to 
be derived from the ISOCAM measurements. Grain tem- 
peratures for the case T c = 3.5 x 10 7 K are also given. 
As seen in the figure the grain temperature of small sili- 
cate and graphite grains is nearly insensitive to T c since 
the electrons are not stopped by the grains. The equi- 
librium temperatures are consistent with the colour tem- 
peratures T c we derived in paper I for a modified Planck 
spectrum Fx oc \~@Bx{T c ), where B\(T C ) is the Planck 
function, ranging from ~ 200 K {ft = 2) to ~ 290 K 
((3 — 0). Because of stronger heating due to the higher 
grain density and the weaker cooling due to lower emissivi- 
ties in comparison to silicate and graphite grains pure iron 
grains attain the highest temperatures. We also note that 
for the conditions in the shocked CSM, graphite grains 
are hotter than silicate grains. In less dense plasmas (e.g. 
njj = 1 cm -3 ) one would expect a slightly higher temper- 
ature for silicate grains (Dwek, 1987, Fischera, 2000). 




0.01 0.10 1.00 

o[/im] 



Fig. 3. Equilibrium temperatures of spherical silicate, 
graphite and iron grains as a function of grain radius a in 
the shocked gas downstream of the blast wave. The grain 
temperatures in model I and II are represented by the 
solid and the long dashed line, respectively. Also shown 
is the effect of raising T from 2 x 10 7 K to 3.5 x 10 7 K 
(short dashed and dotted lines, respectively, for models I 
and II). 

The temperature of small spherical iron grains depends 
strongly on grain size. This was also noted by Chlewicki 
& Laureijs l|1988(l for iron grains heated by the interstellar 
radiation field. 

2.2.2. Grain size distribution in the CSM 

In order to construct a SED for comparison with the data 
a functional form for the grain size distribution must be 
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adopted. The size distribution of the grains in the shocked 
circumstellar environment of SN 1987A is not known and 
might be modified in comparison to the initial distribu- 
tion in the CSM by different processes like evaporation 
during the UV-flash or sputtering in the shocked gas as 
will be shown later. However, for simplicity we assume, 
that the grains in the shocked CSM have a grain size dis- 
tribution similar to grains in the ISM of our galaxy and 
can be described by a simple power law dn oc a~ k da with 
power index k and a minimum and a maximum grain size 
a m i n and a max . Following Biermann & Harwit (19SU|) a 
power law distribution should be a general description of 
a grain size distribution resulting from grain-grain colli- 
sions and would especially describe the grain size distri- 
bution in the atmospheres of red-giants. In contrast to 
the emission from grains in the ISM, heated by the ISRF, 
where only stochastically heated small grains can achieve 
high enough temperatures to emit at shorter wavelengths 
(Draine & Anderson. HMfy . in the shocked CSM of SN 
1987A larger grains can also contribute to the measured 
fluxes. The composite SEDs are relatively insensitive to 
amin, which we fixed at 10 A. 

2.2.3. Free dust model parameters 

For each of model I and model II for the shocked gas 
we made calculations taking various combinations of the 
following parameters as free variables: 

— Md, the dust mass. This was a scaling parameter in all 
calculations. 

— k, the exponent in the grain size distribution dn(a) oc 
a~ k da. 

— a max , the maximum grain size. 

— Grain composition. Five compositions were considered: 
Pure silicate, pure graphite, pure iron, a silicate-iron 
mixture and a silicate-graphite mixture. The relative 
composition in the mixtures was considered as a free 
variable. 

The most probable values for the variables were found 
through a x 2 -fit to the measured flux densities, taken from 
paper I, which we colour corrected (see Blommaert et al., 
2001J on the basis of the modelled spectrum. In addition 
we derived for each fit the luminosity Ld of the theoreti- 
cal dust emission spectrum. The one sigma uncertainties 
of the parameters were calculated from the x 2 -fit by vary- 
ing A xi = Xi - X m in until Ax? = 1 (see e.g. Press et al., 
1992J. To estimate the uncertainties of the fitted parame- 
ters a m i n , k and the dust mixture the dust mass Md was 
taken to be a free variable. For simplicity the uncertainty 
of Md itself was derived with all other parameters fixed. 

3. Results 

The results of the model calculations are summarised in 
Table where we give for each fit also the reduced xi > 
where v is the number of free parameters. The best fits 
obtained for the two mixtures are shown in Fig. ^together 



with the fitted spectrum for silicate grains. The emission 
spectrum is such that measurements at longer wavelengths 
(e.g. Lundqvist et al.. I1999]l are too insensitive to probe 
the emission from this region. The fits can be summarised 
as follows: 

1. Fit parameter: Md. 

We simply compared the measured flux densities in the 
IR with the theoretical emission spectra of three differ- 
ent grain compositions (pure silicate, graphite or iron) 
with fixed grain size distribution for a max = 0.25 /mi 
and k — 3.5 as proposed for the grains in the ISM 

For both model I and model II a better fit is achieved 
under the assumption of pure graphite grains instead 
of silicate grains in the circumstellar environment. As 
seen in Fig.0]the flux at 6.7 /mi is too high to be emit- 
ted from silicate grains. Because of the high predicted 
temperatures a pure iron composition is even more un- 
likely. 

2. Fit parameter: k, Md. 

For each of pure silicate, graphite and iron we var- 
ied the power k of the size distribution with fixed 
a ma x = 0.25 /mi. 

The derived grain size distributions for graphite grains 
in model I and for silicate grains in both models 
are consistent with k = 3.5. For graphite grains in 
model II, which reach too high temperatures, a flatter 
distribution (k < 2.25) is required. No grain size dis- 
tribution could be found which admitted a pure iron 
solution. 

3. Fit parameter: Dust mixture, Md. 

For a fixed grain size distribution with power k = 3.5 
and a max = 0.25 /mi we varied the composition of cir- 
cumstellar grains, for a mixture of silicate with iron 
grains and a mixture of silicate with graphite grains. 
Excellent fits could be achieved for silicate-iron mix- 
ture in model I and for the silicate-graphite mixture in 
model II (see also Fig.0J. 

4. Fit parameter: a max , Mi. 

For a given dust composition (silicate, graphite and 
the two derived mixtures) and a fixed size distribution 
with power k — 3.5 we varied a max . 
The maximum grain sizes in model I are significantly 
smaller than in model II. In model I for silicate, 
graphite or the silicate-graphite mixture a better fit 
can be achieved with a max < 0.25 /mi. The estimated 
maximum grain sizes of all considered compositions 
in model II are several times larger. We also exam- 
ined a pure iron composition and found that the max- 
imum grain size would be unrealistically large with 
a > 10 /mi. 

Due to the large number of free parameters several so- 
lutions with different compositions and grain size distribu- 
tions are possible. The derived values like a max , k or the 
mixture of the composition are interrelated and depend 
on the plasma density. Despite these uncertainties the fol- 
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Fig. 4. Theoretical emission spectra of dust in the shocked circumstellar medium of SN 1987A for model I (hh = 
300 cm~ 3 ) and model II (nH = 600 cm' 3 ) . Shown are spectra from pure silicate (dotted line) and the best fits with 
a mixture of silicate with graphite (dashed line) and with a mixture of silicate with iron (solid line). The grains are 
assumed to have a grain size distribution dn ex a~ 3 5 da with minimum and maximum grain radii of 10 A and 0.25 /im. 
The ISOCAM fluxes (triangles), shown with the widths of the filters as horizontal lines, have been colour corrected 
to the spectrum of the silicate-iron mixture. The shown uncertainties are the absolute uncertainties in the integrated 
flux densities (taken from Table 2 of paper I) which are 25% (6.75 /im), 4.3% (12 /xm) and 5.5% (14.3 /im). These 
uncertainties are a combination of the systematic calibration uncertainties, which dominates at 12 and 14.3 /im, and 
the random uncertainties. 



lowing conclusions could be reached largely independent 
of the assumptions: 

1. Luminosity 

For the luminosity of the collisionally heated dust 11 
years after outburst we obtained (dependent on com- 
position and the assumed plasma density) values for 
in the range 2.91 x 10 28 to 3.35 x 10 28 W. This can be 
compared with the crude estimate of Ld = 2.5 x 10 28 W 
given in paper I where we approximated the dust emis- 
sion spectrum with a simple modified black body spec- 
trum. As already mentioned in paper I the luminosity 
in the MIR is larger than in X-rays. Despite the in- 
crease of the X-ray luminosity since detection ~ 1400 
days after outburst (Hasingcr et al. 1996) the luminos- 
ity in the energy range 0.5 to 10 keV even ~ 13 years 
after outburst was only Lx ~ 2 x 10 28 W (Burrows et 
al.,HIU|). 

2. Composition 

The dust is most likely a mixture of silicate with iron 
or silicate with graphite. Also pure graphite gives a 
reasonable fit. Pure silicate on the other hand seems 
to be unlikely and pure iron can be excluded. 

3. Dust-to-gas ratio 

Considering only cases where a max is not larger than 
0.25 /^m the dust mass is found to be nearly inde- 
pendent of composition and grain size distribution. 
The corresponding acceptable fits for model I (tt-h = 
300 cm~ 3 ) yield dust masses in the range (0.9 to 1.1) x 
10 -6 Mq and those for model II (nn = 600 cm~ 3 ) yield 



Taking these values the dust-to-gas ratio is only mildly 
dependent on gas density as one would also expect. 
The mass of the X-ray emitting gas is given approxi- 
mately by 



rn^EM 



(2) 



dust masses in the range (0.45 to 0.67) x 10 6 M 



©■ 



where EM is the emission measure from X-ray obser- 
vations and ttih the hydrogen mass. We have taken 
the emission measure from Hasinger et al. (1996) with 
(1.4±0.4) x 10 57 cm~ 3 ~ 2500 days after outburst and 
extrapolated this result with t 2m (see appendix |A"|> to 
4000 days. Using equation we derive a gas mass of 
1.0 x 10" 2 /?i H [300 cm' 3 ] M . Overall, this yields a 
dust-to-gas ratio in the range (0.9 to 1.3) x 10 -4 . 
If the grain size distribution is allowed to extend to 
grains with radii much larger than 0.25 /an the ac- 
ceptable fits imply that the dust-to-gas ratio could be 
as high as 2.2 x 10~ 4 . As will be seen in Sect. 0]thc 
pre-supernova dust abundance would however be al- 
most unaffected. 



4. Destruction processes due to the supernova 

To compare the results for dust abundance and dust com- 
position in the CSM of SN 1987A with expectations for 
the CSM prior to the supernova event, account must be 
taken of the grain destruction processes due to the super- 
nova and its remnant. 

Before quantitatively discussing grain evaporation 
(section ^31 and grain sputtering fsection I4.2JI we briefly 
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Table 2. Results from the dust modelling. 




results model I («h = 300 cm a ) 


results model II (nu = 600 cm a ) 


compos. 


mixture 11 k ew[/H M d [lO-°M ] L d [10" B W] X t 


mixture k a m „H MdllO^Mg] L d [10" s W] %l 


single variable: Md (y = 2) 


silicate 

graphite 

iron 


3.5 0.25 0.97 ±0.03" 2.94 6.1 
3.5 0.25 1.04 ± 0.04 3.35 1.6 
3.5 0.25 1.31 ± 0.05 3.35 55. 


3.5 0.25 0.50 ±0.02 2.63 5.3 
3.5 0.25 0.56 ±0.02 3.15 3.5 
3.5 0.25 0.86 ±0.03 3.65 81. 


two variables: Md and k (v = 1) 


silicate 

graphite 

iron 


3.99^0.24 0- 25 0.85 ±0.03 2.83 8.9 
3.56±g:f? 0.25 1.03 ±0.04 3.36 3.1 
< 0.53 0.25 1.51 ±0.05 2.54 7.8 


3.39+.1-.32 0.25 0.51 ±0.02 2.64 11. 

< 2.25 c 0.25 0.67 ±0.02 3.03 1.3 

< 0.28 0.25 1.09 ±0.04 2.59 18. 



two variables: Md and dust mixture (y = 1) 



si.±iron 


2 09 +i ' U9 


3.5 


0.25 


1.11 ±0.04 


3.18 


0.0 


9 QQ + -L-kt> 
z -°°-0.73 


3.5 


0.25 


0.60 ± 0.02 


2.96 


2.4 


si.±iron 


z - au -0.82 


3.5 


0.10 


0.97 ±0.03 


3.18 


0.3 


3.09+^5 


3.5 


0.10 


0.51 ±0.02 


2.95 


4.4 


si.±iron 


+ 1.89 
O.OO_ .95 


3.5 


0.06 


0.93 ±0.03 


3.19 


0.75 


3.39t? : « 


3.5 


0.06 


0.48 ±0.02 


2.96 


5.7 


si.±gra. 


25+ 56 

u - ZJ -0.25 


3.5 


0.25 


1.03 ±0.03 


3.27 


2.4 


0.73+^8 


3.5 


0.25 


0.53 ±0.02 


2.94 


0.2 


si.±gra. 


zn+ - 64 


3.5 


0.10 


0.93 ±0.03 


3.22 


1.07 


o.84tS:2 


3.5 


0.10 


0.47 ±0.02 


2.92 


0.47 


si.±gra. 


u - JZ -0.34 


3.5 


0.06 


0.89 ±0.03 


3.20 


0.55 


0.88± ;I 


3.5 


0.06 


0.45 ± 0.02 


2.91 


1.45 



two variables: Md and a max (y = 1) 



silicate 




3.5 


03 +tJ - U7 


0.78 ± 0.03 


2.78 


9.3 




3.5 


<n+ 1 - m 

U - 01 -0.41 


0.60 + 0.02 


2.70 


10.2 


graphite 




3.5 


iQ+o-30 

u - la -0.15 


1.01 ±0.03 


3.36 


3.1 




3.5 


Q7+1-60 

u - y ' -0.44 


0.75 ±0.03 


3.13 


3.3 


si.±iron 


2.09 


3.5 


0.26tg;g 


1.12 ±0.04 


3.18 


0.0 


2.33 


3.5 


i cc;>10.0 
r.OO_ 1 24 


1.08 ±0.04 


3.05 


0.9 


si.±gra. 


0.25 


3.5 


n 05+ 014 


0.90 ± 0.03 


3.28 


1.0 


0.73 


3.5 


42+ ' 71 


0.59 ±0.02 


2.95 


0.05 



a Relative mass of silicate to iron or silicate to graphite, respectively. 

The uncertainty of Md was derived with all other parameters fixed. 
c The limit given for k corresponds to x 2 W — X 2 (^ = 0) = 1; here Mi, Ld and xt were calculated for k — 0. 



demonstrate that the grains are dynamically coupled to 
the gas due to the betatron effect. For example the po- 
tential of a graphite grain with a = 0.1 fim at a temper- 
ature of Ti = 3 x 10 8 K is more than 30 Volt (Draine & 
Salpeter. I1979a|l . The larmor radius of the grain is given 
by Rl = m gI v gr /Q gr B, where m gI is the mass of the grain, 
v gI the relative velocity of the grain to the gas (initially 
j v s): Qgr the grain charge and B the magnetic field. For B 
we adopt the value for equipartion between magnetic field 
and relativistic particles (Longair,^H|?|), calculated from 
an extrapolation of the synchrotron emission after 1200 
days (Gaensler et al.. H997|) . Assuming the radio emission 
to arise from the same volume as the X-ray emission (see 
appendix |A"|) . this yields B ~ 10 -7 T, in agreement with 
the value found by Ball & Kirk 1)1992(1 . Taking the shock 
velocity to be 2900 km/s, as derived from radio observa- 
tions (Gaensler et al.. H997|l . the larmor radius i?L is then 
of the order of 2.1 x 10 11 m, which is less then ~ 10 -4 of 
the diameter of the radio emission region. By comparison, 
the distance of the outer blast wave to the contact discon- 
tinuity increases with ~ 0.1 vgt (Chevalier, 1X9521 . Thus, 
the grains comoved with the gas after roughly one month. 



To discuss the mass loss due to evaporation and sput- 
tering we assumed, as in our earlier examination of the 
IR emission, that the grains in the CSM at the time 
of the supernova explosion had a grain size distribution 
with k = 3.5. Again we choose as minimum grain size 
Omin = 10 A. 



4.1. Evaporation during the UV-flash 

The evaporation of grains in the neighbourhood of a super- 
nova in general has already been the subject of earlier ex- 
aminations (e.g. Draine & Salpeter. I1979bl Draine. H981I 
Pearce & Mayes, I198(jf) and also has been discussed for 
silicate grains for the SN 1987A (Emmering & Chevalier, 
1989, Timmermann & Larson. Il993f) . All considered the 
effect of evaporation for single grain sizes only. Here we 
are interested in the mass loss of grains with a certain size 
distribution in the shocked plasma behind the outer shock 
front. We use more recent theoretical results for the spec- 
trum, duration and luminosity of the UV flash (Ensman & 
Burrows. I1992|l . Apart from silicate we also consider iron 
and graphite grains. 

As a reasonable assumption, we only consider the evap- 
oration of grains during this UV-flash, when the grains 
reached their highest temperatures. For the luminosity, 
temperature, and duration of the UV-flash we take the 
theoretical results of the model 500fulll from Ensman & 
Burrows ( 1992|). For our calculation, the spectrum is taken 
to be a simple black body spectrum with the given colour 
temperature from Ensman & Burrows normalised to the 
luminosity. 

To simplify the derivation of the mass loss of grains 
during the UV-flash, we neglect the effect of stochastic 
heating and assume that all grains are at their equilibrium 
temperatures. The evaporation of atoms from the surface 
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of a grain with temperature leads to a reduction in 
radius, that can be described by 

da [/H _ 1 10 U e -ff(a)/teT d («) ] / 3 ) 
dt[s] 3 

where fee is the Boltzmann constant and W(a) the energy 
single atoms need to escape from the surface, which de- 
pends on the radius of the grain because of the surface 
tension of the grain. The coefficient is chosen to be close 
to the value given by Guhathakurta & Draine and 
Voit (1991) for graphite and silicate. The energies W(a) 
for silicate and graphite are adopted from Guhathakurta 
& Draine. For iron we adopt 

W{a) = (50 000- 22 000 iV(a)" 1/3 )fc B , (4) 

where N(a) is the number of atoms in the grain. Here 
we have chosen a bounding energy of the atoms close to 
U = 4.29 eV, given in Gerthsen et al. I|1989f) . and a surface 
tension of a = 1.8 J/m 2 , given in Levefre Q979). 

The total reduction a ev (a,t) of a grain with an initial 
radius a during the UV-flash after a time t is found by 
integrating Eq. If the final radius a(a, t) = a — a ev (a, t) 
is smaller than 3 A, the grain is assumed to have evap- 
orated. The final grain size distribution of non evapo- 
rated grains assumed to have an initial size distribution 
dn(a) — f(a)da — A a -3 ' 5 da with constant A given by: 

/(a) da = /(a(a)) -r? da 

= A(a + a ev )- 3 - 5 (l-^) \a. (5) 

The derived size distributions of silicate and graphite 
grains for different distances are shown in Fig. [5j As the 
evaporation is a surface effect, the fractional change in ra- 
dius of grains at the same temperature is larger for smaller 
grains. Therefore the number of small grains decreases 
much more rapidly than the number of bigger grains. In 
addition, at the same distance to the supernova smaller 
grains are generally heated to higher temperatures. Due 
to the exponential dependence of Eq.|31the evaporation of 
smaller grains is faster, which causes the maxima in the 
distribution functions. 

How much silicate, iron and graphite dust might have 
survived the UV-flash in a certain distance to the super- 
nova is shown in Fig.Elfor three different maximum grain 
sizes (0.25, 0.1 and 0.06 /im). It can be seen that silicate 
grains evaporate out to larger distances from the super- 
nova than graphite grains. This is partly due to the higher 
bounding energy of graphite grains but mainly caused by 
the much higher temperatures the silicate grains attain in 
comparison with graphite grains. This is the opposite of 
the situation for collisionally heated dust in the CSM. 

The grains most probably responsible for the measured 
infrared fluxes originated from between the original inner 
boundary of the Hll-region and the position of the blast 
wave 4000 days after outburst. The mass loss of iron and 
silicate grains in this region is comparable and significant 




r/10 ,5 m 

Fig. 6. Survived grain mass of graphite, silicate and iron 
grains after the UV-flash versus distance r from the su- 
pernova. It is assumed, that the grain size distribution 
before the supernova outburst had k — 3.5. a mnl is cho- 
sen to have been 10 A. The three curves (solid, dashed 
and dotted line) for each dust composition correspond to 
different choices for the initial maximum grain size a max 
(0.25, 0.1 and 0.06 /mi). Also shown is the approximate 
range of positions covered by the thick inner ring (dark 
grey) and the Hll-region (light grey) before being caught 
by the blast wave. The inner radius of the Hll-region was 
taken to be 4.72 x 10 15 m. This was calculated from the 
observed position of the blast wave at epoch 3200 days 
(Gaensler et al.,|1997| extrapolated to the epoch at which 
the shock first reached the Hll-region (1200 days after out- 
burst; Staveley-Smith et al.,[1992). The vertical straight 
line at a radius of 5.43 x 10 15 m indicates the position 
of the blast wave at the epoch of the ISOCAM observa- 
tions (4000 days after outburst, also extrapolated from 
the observed position at 3200 days). The extrapolations 
were made assuming a velocity of 2900 km/s for the blast- 
wave corresponding to model II. If we take the velocity of 
4100 km/s (model I) the inner boundary is closer to the 
supernova and the position of the blast wave further out. 



even for the largest considered grain size of 0.25 /im. This 
changes at larger distances from the supernova, where 
evaporation of predominantly small iron grains becomes 
much stronger. For graphite grains evaporation is unim- 
portant for the whole Hll-region. 

The radially integrated evaporated dust masses of sili- 
cate, graphite and iron grains corresponding to the shown 
curves is given in Table |3J In the integration it is assumed 
that the shock surface area is proportional to the square 
of the distance r of the blast wave to the position of the 
supernova. The differences in the derived values for the 
two models are due to the different shock velocities, that 
give slightly different positions of the inner boundary and 
the outer shock. 
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0.001 0.010 0.100 1.000 0.001 0.010 0.100 1.000 
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Fig. 5. Theoretical grain size distributions (dashed lines) of survived silicate and graphite grains after the UV-flash 
at various distances D (given in units of 10 15 m) to the supernova. The initial size distribution is assumed to have 
been a power law dn(a) = v4a" 3 ' 5 da with constant A (straight solid line). 



4.2. Sputtering in the shocked gas 

Behind the blast wave grains undergo sputtering. This is 
thought to be one of the most important destruction pro- 
cesses in fast moving shocks (see e.g. Dwek et al., 1996). 
Sputtering time scales appropriate for a plasma with the 
abundances of the shocked CSM of SN 1987A are given in 
appendix IbI 

The final grain size distribution after a time At of an 
initial grain size distribution /(a) da due to sputtering in 
a hot plasma is given by: 

/(a, At) da = /(a(a, At)) ^ da = /(a + Aa(Ai)) da. (6) 

da 

Because of the evaporation of grains during the UV-flash, 
the grain size distribution f(d) da entering the sputtering 
zone is dependent on the distance to the supernova. The 
average grain size distribution at the time t after the shock 
reached the HH-region at time to is therefore proportional 
to: 

(/(S, t - t )) oc jf dt' v s /(a, t - t', r(t')) r(t'f , (7) 

where to = 1200 days and the surface of the outer shock 
front is again assumed to increase with the square of the 
distance r(t) of the blast wave to the supernova position. 

The derived mass losses for the considered dust com- 
positions and the three different maximum grain sizes are 
tabulated in Tabled! For comparison, we derived the sput- 
tered dust mass with and without previous evaporation. 
The mass loss due to sputtering becomes progressively 
more important for the bigger grains. Whereas sputtering 
reduces the radius independent of the grain size, evapora- 
tion is the dominant destruction process for small grains. 
Graphite grains which are relatively stable against sput- 
tering (appendix[BJ) and did not evaporate during the UV- 
flash (Sect. m suffer only moderate depletions compared 
to silicate and iron grains. In model I the initial mass of 
graphite grains would be less than a factor of 1.3 higher 



Table 3. Mass loss through evaporation and sputtering 



n H a 


300 cm" 3 


600 cm" 3 


£Jmax[M m ] 


0.25 0.10 0.06 


0.25 0.10 0.06 


AMsputt." sil. 
AM cvap . sil. 


27.7% 41.1% 50.3% 
49.4% 80.4% 98.3% 


44.0% 61.7% 71.9% 
48.7% 79.2% 98.9% 


AM S putt. iron 
AM cvap . iron 


27.2% 40.4% 49.5% 
48.0% 78.6% 99.8% 


40.9% 58.1% 68.4% 
47.7% 78.1% 100.% 


AM tot al C Sil. 

AMtotai iron 
AMtotai gra. d 


53.1% 82.4% 98.4% 
51.8% 81.0% 99.8% 
11.3% 17.3% 22.0% 


58.4% 85.1% 99.1% 
56.4% 83.6% 100% 
19.1% 29.0% 36.3% 



" Density of the shocked gas downstream of the blast wave. 
b Sputtering without previous evaporation. 
Mass loss after evaporation and following sputtering. 
For graphite only the total mass loss is given because of an 
insignificant mass loss due to evaporation. 



than that infered from the ISOCAM observations for all 
considered maximum grain sizes. 

5. Discussion 

5.1. Dust in the pre supernvova CSM 

The modelling of the ISOCAM measurements (Sect. |3J), 
combined with the calculations of grain destruction due 
to evaporation and sputtering (Sect. allow us to infer 
grain abundances in the pre supernova CSM as well as to 
make crude estimates of maximum grain sizes and grain 
composition. 

This information was derived for the best solutions 
(xl < 3, Table |2J) using the mass loss estimates given in 
Table 13 and is tabulated in Tabled As expected the main 
parameter determining the composition and abundance of 
grains in the pre-supernova CSM is the maximum grain 
size. It is interesting to compare these quantities with ob- 
served and derived properties of circumstellar dust in stel- 
lar outflows. 
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Table 4. Solutions for dust in the pre-supernova CSM 



model 


comp. 


mixture' 




10" 


4 M d /M gas 






0.25 6 


0.10 


0.06 


0.25 


0.10 0.06 


I 


gra. 








1.1 


1.2 1.3 


I 


si.+iron 


2.15 


3.13 


0.44 


2.3 


5.4 152 


II 


si.+iron 


2.45 






2.8 




I 


si.+gra. 


0.47 


2.02 


25.4 


1.4 


2.4 20. 


II 


si.+gra. 


1.42 


4.01 


62.3 


1.8 


3.6 48. 



a Relative mass of silicate to iron or silicate to graphite. 
b Maximum grain size in microns. 



We are not aware of direct measurements in the IR 
and submm regimes of dust abundances in the winds of 
LMC stars. We will simply estimate it by assuming that 
the dust abundance in the winds is proportional to the 
metallicity of the ISM. This is supported by Woods et al. 
( 1992) who found that the outflow velocity of oxygen rich 
stars in the LMC is significantly lower than for oxygen rich 
stars in our galaxy and suggested that the lower velocity 
is due to the lower dust-to-gas ratio caused by the lower 
metallicity of the LMC. The same has also been proposed 
for high luminosity stars in the galactic anticentre which 
show only a modest outflow velocity (Habing et al.,^94). 
These works were confirmed by van Loon (2000), who 
found by comparing obscured Asymptotic Giant Branch 
stars of our galaxy, the Large and the Small Magellanic 
Clouds that the inferred dust-to-gas ratio of both carbon 
and oxygen rich stars is approximately proportional to the 
initial metallicity. 

The winds of evolved carbon rich stars in our galaxy 
have generally been found to have dust-to-gas ratios in 
the range ~ 0.1 to ~ 1% (Jura, EHHUI Martin & Rogers, 
IT3S7I Griffin, HM1 Knapp et al., Bagnulo et al., 

1995 and Olofsson et al.. 119931 as quoted by Hiriart & 
Kwan l2"nnnjl . On the other hand Hiriart & Kwan pWijl 
estimated from a subsample of Olofsson et al. a maximum 
dust-to-gas ratio of only ~ 0.1%. For the galactic oxygen 
rich star OH 231.8+4.2, Knapp et al. fiM^i obtained a 
dust abundance of ^0.7%. 

Scaling by the factor of 2 between the metallicity of 
the LMC and the metallicity of the solar vicinity (Russell 
& Dopita, 1992) we would expect on this basis dust to gas 
ratios for LMC carbon stars in the range ~ 0.05 — 0.5 % 
and for LMC Oxygen stars of order ~0.3%. Although 
clearly very uncertain, these numbers are nevertheless 
only consistent with a subset of the solutions of the pre- 
supernova gas-to-dust ratios in the CSM of SN 1987A 
given in Table 0] In particular, all the pure graphite solu- 
tions appear underabundant in dust by an order of magni- 
tude compared with expectations for the winds of carbon 
stars. This confirms the expectation that the dust should 
in fact be silicate rich on the basis of the gas phase abun- 
dances in the inner ring (Lundqvist & Fransson, 1996 
Sonneborn et al.. H997jl . 

The solutions for silicate-iron and silicate-graphite 
mixtures in Tableware consistent with the expected dust 



to gas ratios of ~ 0.3% for oxygen-rich stars provided the 
maximum grain size in the CSM of SN 1987A was smaller 
than 0.1 /im. This is in accordance with findings, both 
theoretical and observational, that cool stars eject mainly 
small grains, irrespective of composition. An upper limit 
of 0.14 /mi on the maximum sizes of grains around oxy- 
gen rich mass- losing stars was found by (Jura. H~996|l . A 
maximum grain smaller than ^0.1 /im is also consistent 
with the sizes found for the dust ejected from the car- 
bon rich star IRC+10216 (Martin & Rogers,[THE3 Griffin, 
1990] Jura. 119941 Bagnulo et al.. H995|) or the typical sizes 
around evolved carbon stars derived by Hiriart & Kwan 

mm . 

On the other hand, this small maximum grain size con- 
trasts with results for oxygen rich stars derived from their 
ultraviolet extinction. Rogers, Marting & Crabtree Q1983[) 
suggested that the grains around /i Cep should be in the 
range 0.1 /jm to 0.5 /jm. Based on extinction measure- 
ments made for a Sco, Seab & Snow (1989) concluded 
that the grains in the circumstellar environments of cool 
oxygen rich giants should be larger than ~ 0.08 /im and 
possibly enrich the ISM with grains as large as ~ 1 /im. If 
this were also the case for the RSG wind of the progenitor 
of SN 1987A, the low dust abundance found by ISOCAM 
could neither be explained by evaporation nor by sput- 
tering. It might then have had to have been intrinsically 
low. Alternatively, the grain abundance could have been 
reduced by radiation pressure (Turner & Pearce. ll992[) af- 
ter the progenitor evolved to its final BSG phase 20 000 
years ago (Crotts & Heathcote. I1991|) . However for this, 
the magnetic coupling of the grains to the gas would have 
had to have been weak. Another scenario for reducing the 
grain abundance might be a mixing of the material of the 
BSG wind, with a negligible dust abundance, with that of 
the RSG wind at a dynamically unstable interface between 
the winds (Garcia et al.. I1996al I"l996bj> . 

5.2. Iron abundance in the shocked CSM 

On the basis of a spherically symmetrical hydrodynam- 
ical simulation of the interaction of the blast wave with 
the HII region, Borkowski et al. (|1997f) derived an upper 
limit on the gas phase iron abundance in the HII region 
of only 0.1 of the solar iron abundance from the X-ray 
spectrum measured with the ROSAT satellite (Hasinger 
et al.. H99b1) . Comparing Chandra data with a plane par- 
allel shock model Park et al. ( 2002 ) found a value of 0.07 of 
solar for the iron abundance in the X-ray emitting region 
between epochs 1999 and 2001. These gas phase abun- 
dances for iron are lower by a factor of at least ~ 3 com- 
pared to the iron abundance in the LMC from Russell & 
Dopita (1992), which is 36% of the solar abundance (de- 
termined from the photosphere; Anders& Grevesse 1989). 
This prompted Borkowski et al. to suppose that most of 
the iron was condensed into grains. However, this is not 
supported by the ISOCAM measurements. An upper limit 
for the mass of iron condensed in grains can be taken 
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from the calculations for pure iron grains in the shocked 
CSM tabulated in Table fj] Taking for iron grain masses of 
1.5 x 1O~ 6 M (model I) or 1.1 x 1O~ 6 M (model II) and 
the iron abundance in the LMC to be npe/^H = 1-7 x 10~ 5 
(Russell & Dopita,l]2T32J), the fraction of iron in solid form 
is at maximum ~ 32% (model I) or ~ 45% (model II). For 
the better fitting silicate/iron mixtures the upper limits 
will be still lower. On the basis of this evidence, the HII- 
region would underabundant in iron, whether in gaseous 
or condensed form. Further X-ray and infrared observa- 
tions would be valuable to investigate this problem. 

5.3. Outlook 

In the short term, further observations to follow the MIR 
light curve will provide information on the dependence 
on radial position of the dust-to-gas ratios in the shocked 
HII region for comparison with model predictions. On the 
one hand there will be a tendency for the overall volume- 
averaged dust-to-gas ratios to be lowered due to sputtering 
at later epochs, if the upstream dust abundance is con- 
stant. On the other hand, the survived dust abundance 
after the UV- flash will increase with radius, especially for 
silicate and iron grains. Another potential reason for in- 
creasing dust abundances with radius could be a flushing 
out of grains from the inner regions of the CSM through 
radiation pressure, after the RSG turned into a BSC This 
might also offer an explanation for the puzzle of the low 
iron abundance in both grains and gas inferred from the 
ISOCAM and ROSAT and CHANDRA data, as discussed 
above. As already stated this may require a weak magnetic 
coupling of the grains to the gas. 

Observations of the brightening so-called "hot spots" 
seen in the optical (Lawrence et al.. I2000|l show that the 
blast wave is already interacting at certain places with 
dense material of the thick inner ring. The origin of the 
thick inner ring, is unknown. It is also not clear whether 
the MIR echo seen after 580 days (Roche et al., I1993[l 
can be attributed to the ring. If the ring is composed of 
material from the red supergiant phase of the progenitor 
star (see e.g. Fransson et al.. 11989)) . then one might antic- 
ipate a rapid increase in MIR luminosity with time, which 
would soon dominate the continuum emission from the 
ffll-region. The higher gas densities in the ring might lead 
to an accompanying increase in grain temperature, which 
would allow thermal dust emission from the thick inner 
ring to be distinguished from an increasing contribution 
from the HII region. If, on the other hand, the thick in- 
ner ring is composed of material from the companion star 
in the putative binary system (e.g. Podsiadlowski. IT992f) 
then one might speculate that the dust abundance and 
composition of the ring might deviate markedly from that 
of the HII region discussed in this paper. 

6. Summary 

The IR emission was analysed in terms of thermal emis- 
sion from dust, collisionally heated by the shocked gas 



behind the blast wave that is expanding into the HII- 
region interior to the thick inner ring. We used a realistic 
grain model with a grain size distribution dn oc a~ k da 
including stochastically heating of small grains. For the 
shocked gas we considered two different models corre- 
sponding to the downstream densities nn = 300 cm -3 
(model I) and tin = 600 cm" 3 (model II). The conclu- 
sions for the shocked CSM are as follows: 

1. The luminosity of the dust emission is found to range 
from 2.91 to 3.35 x 10 28 W (according to assumed grain 
composition and plasma density). 

2. The dust is most likely composed either of a silicate- 
iron or a silicate-graphite mixture or of pure graphite. 

3. Considering only cases with a maximum grain size not 
larger than 0.25 /mi the dust masses are, independent 
of the assumed grain composition and size distribution, 
(0.9 to 1.1) x 10" 6 M (model I) and (0.45 to 0.67) x 
10 -6 (model II). This corresponds for a gas mass of 
1.0 x 10~ 2 /n H [300 cm~ 3 ] M to dust-to-gas ratios in 
the range (0.9 to 1.3) x 10 -4 . For larger maximum 
grain sizes a larger range of the dust-to-gas ratio of up 
to 2.2 x 10~ 4 are admitted by the data. 

4. For the LMC abundances of Russell & Dopita ifTMS)) 
the maximum fraction of iron condensed into grains in 
the shocked CSM is 32% (model I) or 45% (model II). 

We have calculated the grain destruction due to the UV- 
flash of the supernova outburst and subsequent sputtering 
downstream of the blast wave. From this we have deduced 
dust properties in the CSM prior to the supernova. Our 
conclusions are: 

1. Pure graphite solutions for the pre-supernova CSM 
are underabundant in dust by an order of magnitude 
compared with estimates for the abundance of carbon 
grains in the winds of LMC-carbon stars. 

2. Solutions for silicate- iron and silicate-graphite mix- 
tures in the pre-supernova CSM are consistent with 
the expected dust-to-gas ratios of ~ 0.3% in the winds 
of oxygen rich LMC stars provided the maximum grain 
size was less than ~ 0.1 /xm. 
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Appendix A: X-ray flux until 4000 days after 
outburst 

Measurement of the X-ray flux from the SN 1987A un- 
til ~ 3000 days after outburst with the X-ray satel- 
lite ROSAT (Hasinger et al.. I1996|l showed a monotonic 
increase of the luminosity, where the variation in time 
could be linear or also steeper with cx t 1 - 67±0 ' 35 . A lin- 
ear trend would be consistent with a constant external 
density (Hasinger et al.. H996[> . Later measurements with 
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CHANDRA ~ 13 years after outburst gave fluxes clearly 
above a linear trend (Burrows et al., 2000). To see how 
the X-ray flux evolved until the 1SOCAM observations 
we derived the X-ray fluxes from a large number of avail- 
able measurements that where made with the HRI and 
the PSPC instrument of ROSAT until - 4000 days after 
outburst. For a direct comparison with the data published 
by Hasinger et al. we also included measurements which 
were made before ~ 3000 days after outburst. The data 
were automatically calibrated with SASS. In calculating 
the counts of X-ray photons we used the same apertures to 
determine the source and the background counts used by 
Hasinger et al. (1996 ) and scaled the photon fluxes derived 
from the HRI-data with a factor 2.65 to allow comparison 
with the PSPC-data. To estimate the uncertainties we as- 
sumed poisson statistics for the source counts and gaussian 
statistics for the background noise. 



Table A.l. Derived rates of X-ray photons towards 
SN 1987A. 



obs. 


time 


day 


instr." 


iint/s 


N/1000 s 


12.02.91 


- 13.02.91 


1448 


HRI (1) 


23107 


0.17 ± 0.90 


06.10.91 


- 07.10.91 


1685 


PSPCB 


16398 


2.72 ± 0.63 


30.04.92 


- 14.05.92 


1898 


PSPCB 


9340 


2.05 ± 0.76 


04.12.92 


- 06.12.92 


2110 


PSPCB 


2552 


5.12 ± 1.80 


07.04.93 


- 10.04.93 


2235 


PSPCB 


11259 


2.63 ± 0.73 


20.06.93 


- 05.07.93 


2315 


PSPCB 


10391 


3.13 ± 0.76 


28.09.93 


- 30.09.93 


2409 


PSPCB 


9131 


3.79 ± 0.85 


20.06.94 


- 20.09.94 


2718 


HRI (4) 


12223 


5.01 ± 1.62 


03.10.94 


- 02.01.95 


2823 


HRI (2) 


18756 


3.54 ± 0.91 


01.04.95 


- 11.07.95 


3008 


HRI (3) 


28833 


4.61 ± 0.97 


10.10.95 


- 10.01.96 


3196 


HRI (2) 


26034 


7.29 ± 1.15 


16.04.96 


- 31.07.96 


3392 


HRI (2) 


46907 


7.54 ± 0.89 


22.10.96 


- 12.01.97 


3568 


HRI (2) 


45192 


7.22 ± 0.87 


21.02.97 


- 02.03.97 


3654 


PSPCB 


15440 


8.31 ± 0.86 


04.03.97 


- 03.06.97 


3706 


HRI (2) 


52601 


8.31 ± 0.86 


16.12.97 


- 17.12.97 


3948 


HRI (1) 


21738 


12.6 ± 1.5 


19.02.98 


- 22.02.98 


4014 


PSPCB 


19382 


9.19 ± 0.78 



a Number in brackets give the number of observations used. 



The derived final counting rates of the observations are 
summarised in table IA.1I and shown in Fig. IA.1I as open 
symbols. Until 3000 days after outburst they are consis- 
tent with the counting rates published by Hasinger et al. 
which are shown as black solid symbols. At later times 
the fluxes lay slightly above the linear trend. Fitting a 
potential to the counting rates derived here gives with a 
Xmin = 1-29 a monotonic increase of: 

N(t) oc t 2m±Q - 20 . (A.l) 

This may be because the shock is reaching denser material 
in time. On the other hand, numerical calculations for a 
HH-region with homogeneous density do not give a linear 
trend of the X-ray luminosity for the first several years 
(Borkowski et al.. H997|l . 

To derive the volume V of the X-ray emitting gas at 
the time of the ISOCAM observations ~ 4000 days after 
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Fig. A.l. Time evolution of the X-ray emission of SN 
1987A, measured with the ROSAT satellite. The pho- 
ton fluxes derived from observations with the HRI and 
the PSPC instrument are shown as triangles and squares. 
Published data (filled symbols) until 3000 days after out- 
burst (Hasinger et al., I1996f) are shown for compari- 
son. These data can be described through a linear in- 
crease (dotted line) since 900 days and a potential (bro- 
ken line) with oc < 167 (Hasinger et al., 1996). The photon 
fluxes derived here (open symbols) give an increase with 

OC <2.06±0.20 (solid Hne) 

outburst we assumed that the X-ray luminosity was pro- 
portional to the emission measure EM — n c riiV, where 
n e and n\ are the number density of the electrons and the 
ions. Taking the emission measure after ~ 2500 days to be 
(1.4 ± 0.4) x 10 57 cm -3 (Hasinger et al.. 119961 Borkowski 
et al.. 11997)) and the abundances as given in Sect. 12. ll this 
volume was approximately 2.2 x 10 52 cm 3 . 

Appendix B: Sputtering time scales for the 
shocked CSM 

Here we derive sputtering rates specific to the physical 
conditions and abundances encountered in the shocked 
HII region around SN 1987A. At the high temperatures 
involved, sputtering is independent of grain charge and 
nearly independent of the plasma temperature. Sputtering 
reduces the radius of all grains of a given composition in 
a time At by the same amount Aa(At). To derive the 
life times of the grains due to sputtering we used the for- 
mula (Eq. 27) given by Draine & Salpeter (1979a) for non 
rotating grains moving with a relative velocity v gr with re- 
spect to the gas with a temperature T{. Sputtering yields 
were adopted from Tielens et al. I|1994|l . and we consid- 
ered sputtering due to H, He, C, N and O, adopting plasma 
abundances as given in Sect. 12. ll The initial relative veloc- 
ities of grains overtaken by the blast wave (v gT = It's) are 
so high that at first the sputtering is almost non thermal. 
Due to the drag forces the sputtering becomes thermal af- 
ter some time, which increases the sputtering yield slightly 
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(10%-30%). Here we simply take the average of the sput- 
tering yield of the two limits v gI = and v gI — 3u s /4. The 
resulting life times of the considered grains arc: 



_ o[0.0lA*m] 



in model I and 



_ o[0.01/im] 

)iH [cm -3 ] 



in model II. 



11.5 x 10 3 graphite 
2.3 x 10 3 silicate years (B.l) 
2.7 x 10 3 iron 



12.8 x 10 3 graphite 
2.9 x 10 3 silicate years (B.2) 



3.0 x 10 3 



iron 
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